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4.  Description  of  the  Research  Progress 

4.1  Ferroelectric  thin  film  synthesis  and  texture  analysis 

Sol-gel  deposition  and  RF  sputtering  process  was  developed  for  deposition  of  PZT  on  Pt/Ti/Si02/Si  (hereafter, 
referred  to  as  platinized  Si)  substrates  as  shown  in  Fig.  1.  Target  preparation  was  perfected  for  Zr/Ti  ratios  of 
60/40  and  52/48  with  excess  Pb  to  compensate  for  Pb  loss  during  post  deposition  annealing.  As  deposited  PZT 
RF  (henceforth,  RF  refers  to  “RF  sputtering”  in  this  report)  thin  films  were  not  well  textured  (i.e.  with 
preferred  crystalline  orientation).  To  texture  and  obtain  crack-free  thick  PZT  RF  films,  we  employed  pre¬ 
treated  substrates  and  post-deposition  annealing.  One  pre-treatment  was  the  use  of  seed  layer  of  textured  PZT 
sol-gel  thin  film  of  thickness  65-85nm  [1]. 


Figure  1:  Sol-gel  process  flow  and  the  RF  sputter  configuration  for  deposition  of  PZT  thin  films. 

A  detailed  study  was  conducted  to  determine  the  conditions  for  obtaining  preferred  crystallite  orientations 
(referred  as  textured).  The  results  of  this  study  were  summarized  in  Temperature-Time-Transformation  (TTT) 
diagrams.  These  diagrams  provide  two-dimensional  relationships  of  crystalline  orientation  to  pyrolysis  and 
annealing  conditions.  To  augment  our  understanding  of  the  thermal  budgets  required  for  the  texturing  of  the 
PZT  sol-gel  thin  films,  we  further  developed  relationships  between  each  phase  and  the  experimental 
conditions.  In  addition,  the  optical  band  gap,  morphology  and  composition  of  highly  textured  sol-gel  thin 
films  were  evaluated  using  Variable  Angle  Spectroscopic  Ellipsometry  (VASE),  Raman  scattering, 
piezoresponse  force  microscopy  (PFM)  and  X-ray  photoelectron  spectroscopy  (XPS).  RF  sputter  deposition 
operating  space  was  explored  using  statistically  designed  experiments  (using  Design  Expert  software)  and 
ANOVA  (Analysis  of  Variance)  models  were  formulated.  The  responses  were  thickness,  refractive  index  and 

absorption  coefficient  from  ellipsometric  data  and  the  elemental  compositions  from  energy  dispersive  X-ray 
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analysis  of  the  PZT  thin  films.  These  models  were  then  utilized  to  prediet  and  optimize  the  proeess  eonditions 
neeessary  to  obtain  the  preferred  responses. 


Temperature-Time-Transformation  model  for  sol-^el  deposited  PZT  (60/40)  films:  We  eondueted  detailed 
investigation  on  the  texture  evolution  as  this  will  direetly  influenee  the  eleetrieal  response  of  the  high 
frequeney  sensor  strueture.  Based  on  our  statistically  designed  experiments  the  two  parameters,  annealing 
temperature  and  time,  were  more  significant  than  the  pyrolysis  conditions  in  achieving  proper  texture.  This  is 
clearly  evident  in  the  half  normal  probability  plots  of  the  Effects  (measure  of  the  process  variable’s  influence 
on  the  response)  shown  in  Fig.  2.  The  responses  were  normalized  XRD  peak  heights  for  [100],  [110]  and  [111] 
orientations  at  2-theta  angles  of  22°,  31°  and  39°  respectively. 
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Figure  2:  Half  normal  probability  plots  of  XRD  responses  -  [100],  [110]  and  [111]  peak  heights. 


Based  on  this  information,  further  extensive  exploration  of  the  sol-gel  thermal  budget  operating 
space  was  conducted.  These  experiments  are  summarized  in  Temperature-Time-Transformation  (TTT) 
diagrams  shown  in  Fig.  3.  This  is  similar  to  the  diagram  that  was  reported  by  Chen  and  Chen  [2]  except  that 
their  main  process  variables  were  pyrolysis  temperature  and  time  and  they  had  maintained  constant  annealing 
temperature  and  time.  However,  we  wanted  to  expand  beyond  these  pictorial  guides  and  decided  to  develop 
mathematical  relationships  between  XRD  peak  data  and  thermal  processing  conditions.  The  data  utilized  to 
create  the  TTT  diagrams  were  analyzed  using  JMP  statistical  software  and  the  quadratic  fits  are  shown  in  Fig. 
4.  Despite  moderate  R^  values,  the  predictability  of  the  quadratic  model  was  found  to  be  poor.  A  good  hint 
towards  this  unpredictability  can  be  witnessed  in  the  portion  of  points  that  do  not  trend  with  the  surface  plots. 
The  main  reason  for  this  poor  fit  was  that  we  did  not  take  into  account  proportionality  between  different 
crystalline  orientations  (see  Scatter  plot  matrix  correlation  in  Figure  5). 


NoPyrolyjls 


PYrDlysis-25iD  IkS  minutes 


□  (IDOi 


Figure  3:  Temperature-Time-Transformation  diagrams  of  PZT  sol-gel  thin  films  pyrolyzed  at  a)  No  pyrolysis, 
b)  250°C,  1.5  minutes,  and  c)  300°C,  3  minutes. 
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Figure  4:  JMP  contour  plots  of  PZT  sol-gel  thin  films  pyrolyzed  at:  a)  none,  b)  250^C  for  L5min,  and 
300^ C  for  3min,  prior  to  annealing. 
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100  110 

Figure  5:  Scatter  plot  matrix  of  the  responses  (XRD  peak  data). 

We  next  utilized  Logistic  regression  of  the  XRD  peak  data  against  the  3  process  variables  -  pyrolysis 
temperature/time,  annealing  temperature  and  annealing  time.  We  lumped  the  pyrolysis  variables  of 
temperature  and  time  together  as  they  were  blocked  experiments  with  both  variant  between  the  blocks.  We 
therefore  in  statistical  terms  made  pyrolysis  temperature/time  as  categorical  variables.  The  normalized  XRD 
peak  data  response  for  three  different  crystalline  orientations  of  [100],  [110]  and  [111]  were  easily  converted 
to  proportions.  The  final  statistical  responses  were  probabilities  of  occurrence  of  each  orientation  in  a 
particular  thin  film  sample  at  the  stipulated  thermal  conditions.  These  probabilities  are  legit  transformation  of 
the  unordered  proportions  of  [100],  [110]  and  [111].  These  responses  were  modeled  against  a  linear  predictor 
model  of  the  thermal  categorical  and  continuous  factors.  Unlike  pyrolysis  variables  which  were  fixed  as 
categorical,  we  attempted  both  continuous  and  categorical  versions  of  the  annealing  temperature  and  time 
variables. 

Probability [100]  =  1/(1  +  Exp(  -Lin[100]  )  +  Exp(Lin[110]  -  Lin[100])) 

Probability [110]  =  1/(1+  Exp(Lin[100]  -  Lin[110] )+  Exp(-Lin[110])) 

Probability  [111]  =  1/(1+  Exp(  Lin[100])+  Exp(Lin[110])) 

=  1-  Probability [110]  -  Probability [100] 

For  categorical  predictor  factors, 

Lin[xxx]  =  Intercept  +  function  of  (pyrolysis  time,  annealing  temperature,  pyrolysis 

time"^ annealing  temperature,  annealing  time,  pyrolysis  time"^ annealing  time) 

For  continuous  predictor  factors, 

Lin[xxx]  =  Intercept  +  function  of  (pyrolysis  time,  annealing  temperature,  pyrolysis 

time"^  annealing  temperature,  annealing  time,  pyrolysis  time"^  annealing  time,  annealing 
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temperature'^annealing  time,  pyrolysis  time"^ annealing  temperature'^annealing  time) 


For  prediction  when  using  categorical  predictor  factors,  only  annealing  factor  values  used  in  model  were 
permitted  whilst  for  continuous  predictor  factors  any  annealing  factor  values  within  range  were  permitted. 
We  created  4  samples  shown  in  Table  1,  SGI,  SG2,  RVl  and  RV2  where  SG  and  RV  refers  to  PZT  sol 
mixed  up  by  two  different  operators  but  each  at  four  different  pyrolysis  and  annealing  conditions. 


Table  1:  XRD  normalized  data  vs.  model  predictions  for  four  different  samples. 


Sample 

Categories 

Actual 

Categorical  Model 

Continuous  Model 

SGI  1 

[100] 

0.88 

0.924 

0.878 

[110] 

0.12 

0.076 

0.0417 

1  [111] 

0 

4.29E-06 

0.08 

SG2 

[100] 

0.896 

0.924 

0.878 

1  [110] 

0.104 

0.076 

0.0417 

[111] 

0 

4.29E-06 

0.08 

RVl 

[100] 

0.91 

0.924 

0.878 

[110] 

0.09 

0.076 

0.0417 

n  [111] 

0 

4.29E-06 

0.08 

RV2 

[100] 

0.042 

0 

0 

1 

1  [110] 

0.083 

0.022 

0.034 

[111] 

I  0.875 

1 

0.978 

0.966 

Mad^l  Vatidaticin 


Figure  6:  Comparison  of prediction  results  for  four  sol-gel  samples  shown  in  Table  1  -  red  is  Actual,  Green  is 
Categorical  and  Purple  is  Continuous  fit. 

We  observed  that  the  data  utilized  to  create  the  Temperature-Time-Transformation  (TTT)  diagrams  do  not 
follow  any  simple  multiple  regression  (like  linear  or  multi  order  polynomial;  Figure  7))  trends.  The  reason 
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being  that  material  textured  in  one  phase  or  orientation  will  have  lesser  proportion  of  the  other  phases  or 
orientation.  For  example,  PZT  52/48  has  three  dominant  textures,  <100>,  <11 0>  and  <11 1>  and  therefore  the 
TXT  data  for  sueh  a  system  is  trinomial  and  interdependent  (Figure  8). 

•  No  order  or 


Anrte*T(l 

Tempemur  Pynif^  TBfTwrBtiijf 

1®  no  til  *  t™  f 


Figure  7:  Scatter  plot  matrix  showing  the  lack  of  a  linear  or  higher  order  trends  between  PZT  crystalline 
orientations  and  TTT  conditions. 


Figure  8:  Scatter  plot  matrix  of  the  responses  (XRD  peak  data). 

We  utilized  two  multivariate  regression  approaehes,  multinomial  logistie  and  log  ratio  multivariate 
regressions  of  the  XRD  peak  data  against  the  3  proeess  variables  -  Pyrolysis  Temperature/Time  and 
Annealing  Temperature  and  Time.  The  eonstitutive  equations  are  given  below. 

Lln[100]  =  In  =  xa 

VP[iii]/ 

Lin[110]  =  =  xp 

P[ioo]  =  ( 1  +  exp(— Lin[100])  +  exp(Lm[110]  —  Lm[100]))”^ 
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P[iio]  =  (  1  +  exp(-Lm[110])  +  exp(Lm[100]  -  Lm[110]))  ^ 


P[iii]  =  (  1  +  exp(Lm[100])  +  exp(  Lm[110]))  ^ 


LPl  =  In 


[100] 

iinj 


LR2  =  In 


[110] 

[111] 


[100]  =  (1  +  exp(-LPl)  +  exp(LP2  -  LRl))-^ 


[110]  =  (1  +  exp(-LP2)  +  exp(LPl  -  LR2))-^ 


[111]  =  (1  +  exp(LPl)  +  exp(LP2))-i 


For  Categorical  Predictor  Factors, 


Xa  =  tto  +  ttpP  +  axmTm  +  ap^xmP  *  Tm  +  ap^xpP  *  Tp  +  axpTp 


For  Continuous  Predictor  Factors, 


XP  —  Po  +  PpP  +  PrinTm  +  Pp*TmP  *  Tm  +  Pp^xpP  *  Tp  +  PxpTp 


Xa  =  tto  +  ttpP  +  axmTm  +  ap=,xmP  *  Tm  +  ap^xpP  *  Tp  +  axpTp  +  axm*TpTm  *  Tp  +  ap=,xm*TpP  *  Tm  *  Tp 
XP  —  Po  +  PpP  +  PrinTni  +  PphcThiP  *  Tm  +  Pp=,xpP  *  Tp  +  PxpTp  +  Pxm*TpTm  *  Tp  +  Pp*xm*TpP  *  Tm  *  Tp 


M^ilthioniidl  Soitui  Mted  Log  R^^tio 
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Figure  9:  Multivariate  regression  fits  to  TTT  data. 

From  Figure  9,  we  observe  that  multinomial  categorical  gave  the  best  fit  to  the  TTT  data.  For  model  validation, 
we  created  4  samples,  SGI,  SG2,  RVl  and  RV2  and  the  model  predictions  from  multinomial  categorical 
approach  was  confirmed  as  the  best  (Figure  10) . 
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Figure  10:  Comparison  of  prediction  results  for  4  sol-gel  samples  shown  in  Table  1  -  red  is  Actual,  Green  is 
Categorical  and  Purple  is  Continuous  fit. 


4.2  Interfacial  and  structural  studies  of  textured  PZT  sol-gel  samples 

Initially  commonly  used  VASE  analysis  was  utilized  to  study  the  PZT  sol-gel  films  on  platinized  Si 
substrates.  But  the  analysis  was  complicated  by  poor  fit  between  the  model  and  data  (mean  square  error  values 
in  Fig.  11).  Film  inhomogeneities  were  accounted  in  the  model  and  better  fits  were  obtained  by  utilizing 
optical  non-idealities  (non-uniformity  in  measurement  spot). 


Figure  11:  VASE  modeling  issues  on  PZT  sol-gel  thin  films:  (a)  as-is  MSE  of  67.02,  (b)  MSE  of  35.42  with 
non-ideal  modeling,  and  (c)  MSE  of  11.36  with  aperture. 


However,  further  investigation  revealed  that  as  VASE  measurement  setup  utilizes  an  optical  beam  incident  at 
an  angle  to  the  surface  (Fig.  12),  we  were  sampling  the  edge  bead  (see  Fig.  12)  of  the  spin  coated  sol-gel  films 
and  excluding  this  non-uniform  film  around  the  sample  edges,  lead  to  better  fits.  The  fix  was  the  combination 
of  lower  angle  of  incidence,  use  of  an  aperture  as  shown  in  Fig.  13  and  higher  data  averaging. 


1 .  Known  input 
polarization 


^  woouam  co.,  inc. 


Figure  12:  Angle  of  incidence  of  assaying  light  in  VASE. 


X 


X 


spot  Size  -  Ilium 
^  65'T5  degree  ange  of 
rncidence 

+  DepolarizMion  from  sampling  of 
edge  non-uniformity 


Focused  spot  size  =  5mm 
!  55-65  degree  ange  of 
incidence 

4  Aperture  -  6mm  (hole  punched 
paper  square) 

1  Averagng  to  improve  S/N  ratio 

Fbcu£»d  spot 


P3per  Aperture 

— ^itli  non-unitormity 


/ 

V 

J 

/ 

Figure  13:  Use  of  an  aperture  resolved  the  lack  of  fit  issues. 

To  understand  the  effect  of  substrate  on  heterogeneous  nucleation  of  PZT  sol-gel  thin  films,  we  conducted 
depth  profiling  using  incremental  Ar  sputtering  followed  by  X-ray  photoelectron  spectroscopy  as  shown  in  Fig. 
14.  We  observed  an  artifact  of  the  XPS  depth  profiles  (knock-on  effect),  that  is  there  was  an  unnatural  large 
gradient  of  elemental  concentrations  at  the  interface  between  the  PZT  film  and  platinum  under  layer.  In  an 
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85nm  PZT  thin  film,  20-30  nm  of  grading  would  be  easily  diseriminated  by  the  VASE  analysis  and  therefore 
this  gradient  was  unreal. 


Figure  14:  XPS  depth  profile  of  PZT  sol-gel  sample. 


We  eondueted  the  AFM  analysis  of  the  sputtered  erater  as  shown  in  Fig.  15  and  found  20-30nm  roughness  that 
obliterates  any  sharp  transition  at  the  PZT  and  Pt  interfaee.  High  resolution  binding  energy  XPS  depth  profile 
seans  shown  in  Fig.  16  were  then  attempted  and  we  were  unable  to  diseem  any  interfaeial  heterogeneity. 


Figure  15:  20nm  roughness  at  PZT-Pt  interface  in  sputtered  crater  of  the  XPS  depth  profiled  sol-gel  sample. 


Figure  16:  High  resolution  binding  energy  XPS  depth  profiling  of  PZT  sol-gel  thin  film. 
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Ellipsometric  data  can  be  used  to  determine  the  optical  band  gap  of  highly  textured  PZT  thin  films  using  a 


Tauc  plot  (Fig.  17)  for  direct  band  gap  material  where  a  = 


4Uk 


is  computed  from  absorption  coefficient  (k) 


and  its  square  is  plotted  against  incident  light  energy  [3].  The  x-axis  intercept  of  the  tangent  to  the  sloped  line 
beyond  the  absorption  edge  gives  us  the  optical  band  gap  Eg.  We  see  the  Eg  trends  as  (111)  >  (100)  >  (110). 
The  triangular  area  between  the  Tauc  curve  and  the  tangent  is  a  measure  of  the  sub-band  gap  density  of  trap 
states  and  this  trend  was  found  to  be:  (110)  >  (111)  >  (100). 


(ahv)^ 


0.014 


Energy  (eV) 


Figure  1 7:  Tauc  s  plot  to  determine  the  optical  band  gap  of  highly  textured  PZT  thin  films. 


The  deposition  process  variables  were  -  argon  and  oxygen  flows,  chamber  pressure,  RE  power  (DC  Bias), 
working  distance,  temperature,  rotation  speed  and  substrate  bias.  A  statistically  designed  experiment 
(SDE/DOE)  for  8  variables  will  consist  of  2^  =  256  runs  (#  of  runs  =  2"^  where  n  is  #  of  variables).  This 
experimentation  is  best  conducted  in  one  session  to  limit  the  variance  and  such  a  large  number  of  runs  will 
require  considerable  sample  preparation  and  time  allocation.  So  it  is  best  advised  to  decrease  the  number  of 
variables  or  combine  them  in  the  first  screening  experiment.  We  combined  O2  and  Ar  into  varying  O2  flow  in  a 
fixed  Ar  flow,  fixed  working  distance  and  rotation  speed  and  blocked  for  substrate  temperature  at  two  settings 
of  room  temperature  and  lOO^C.  In  doing  so,  we  reduced  the  experiment  to  2"^  =  16  runs  +  2  temperature 
blocks  +  2  center  runs  per  block  =  TOTAL  of  18  experiments  (Table  2). 


Figure  18:  Half  normal  probability  plots  of  VASE  responses  -  thickness,  refractive  index  and  absorption 
coefficient. 
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Table  2:  Full  factorial  statistical  design  of  experiment  in  four  process  variables  and  the  measured  VASE 
responses. 


Std 

Block 

Factor  1 

A;  At:  02  ratio 

Factor  2 
a  Pressure 
mT 

Factor  3 
C:RF  Power 
W 

Factor  4 
D:RF1  bias 
W 

Response  1 
thickness 

A 

Response  2 
n@633nm 

Response  3 
k@633nm 

8 

1 

Hig(h  Temperature 

4.00 

7.00 

200.00 

OOO 

761.92 

1  78^ 

0,0038439 

14 

2 

Temperature 

4.00 

3.00 

200.00 

10.00 

399.17 

1.8497 

0.01 3923 

5 

3 

Htgih  TelfiSei'^'^e 

2.00 

3.00 

200,00 

0,00 

1308,32 

22545 

2.052E-005 

9 

4 

Hi^  Temperature 

2.00 

3.00 

100.00 

10.00 

376.88 

2.0173 

0.018048 

15 

5 

Hi^  Temperature 

2.00 

7.00 

2mm 

10.00 

662.38 

1.9565 

1.6201E-010 

3 

6 

Higih  Temperature 

2.00 

7.00 

100.00 

0,00 

668.47 

22757 

0  34823 

17 

7 

Higih  Temperature 

3.00 

5.00 

150.00 

5.00 

393.99 

1.6314 

□  01 2353 

2 

8 

Temperature 

4,00 

3,00 

100,00 

0.00 

471 .82 

1  8707 

0,096903 

12 

9 

High  Temperature 

4,00 

7,00 

100,00 

10,00 

No  deposition 

18 

10 

Room  Temperature 

3.00 

5.00 

150.00 

5.GO 

426.24 

1.8802 

0.021286 

16 

11 

Room  Temperature 

4.00 

7.00 

200.00 

10.00 

442.44 

1.9399 

0.026216 

4 

12 

Room  Temperature 

4.00 

7.00 

100.00 

OOO 

3S4  42 

1,61 24 

0.033182 

6 

13 

Room  Temperature 

4.00 

3.00 

200.00 

000 

951  25 

20639 

0.016531 

11 

14 

Room  Temperature 

2.00 

7.00 

100,00 

10.00 

355.04 

1.7938 

0,0061854 

10 

15 

Room  Temperature 

4,00 

3,00 

100,00 

10.00 

No  deposition 

13 

16 

Room  Temperature 

2.00 

3.00 

200.00 

10.00 

758.56 

2.0214 

1 .4639E^011 

7 

17 

Room  Temperature 

2.00 

7.00 

200.00 

0.00 

1369.34 

2.1445 

0.00062516 

1 

18 

Room  Temperature 

2,00 

3,00 

100,00 

0.00 

564.84 

2.5352 

0.57916 

The  responses  were  VASE  measured  thin  film  parameters  -  thickness,  refractive  index  (n  @  633nm)  and 
absorption  coefficient  (k  @633nm).  Despite  having  two  thin  film  deposition  runs  being  immeasurable,  the 
significant  factors  were  easily  decipherable  as  shown  in  the  Half  Normal  Probability  plots  in  Fig.  18.  The 
results  of  the  analysis  of  variance  of  VASE  data  gave  us  the  following  models  for  responses  (refractive  index 
and  absorption  coefficient  were  transformed  for  more  linearity  in  the  model). 

Thickness  =-158.01581+906.09265*  02  :Ar  ratio+7.51601*  RF  Power+21.59485*  RFl  bias-14.30249* 
02  :Ar  ratio  *  RF  Power+100. 05294*  02  :Ar  ratio  *  RFl  bias-0.43622*  RF  Power  *  RFl  bias 

1.0/(n@633nm)  =  +0.47896+0.21585*  02  :Ar  ratio 

Ln(k@633nm)  =  +11.25209-21.87405*  02  :Ar  ratio-0.75347*  Pressure-0.11680*  RF  Power+0.71603* 

RFl  bias-2.33929*  02  :Ar  ratio  *  Pressure+0. 25293*  02  :Ar  ratio  *  RF  Power- 13. 4573 7*  02  :Ar  ratio  *  RFl 
bias+6.79026E-003*  Pressure  *  RF  Power-0.010910*  RF  Power  *  RFl  bias+0. 098354*  02  :Ar  ratio  *  RF 
Power  *  RFl  bias 

[Pb]  =-5.15252+17.71069  *  02  :Ar  ratio+0.092922  *  RF  Power+0.34365*  RFl  bias-0.28374  *  02  :Ar  ratio  * 
RF  Power+3. 57606*  02  :Ar  ratio  *  RFl  bias-7. 58598E-003*  RF  Power  *  RFl  bias 

Ln(Ti)  =+1.09901-2.93732  *  02  :Ar  ratio 

[Zr]  =-0.12168+5.48020  *  02  :Ar  ratio-0.13771  *  Pressure+0.025982  *  RF  Power+0.12614  *  RFl  bias- 
0.10126  *  02  :Ar  ratio  *  RF  Power+1. 12925  *  02  :Ar  ratio  *  RFl  bias-1.913 14E-003  *  RF  Power  *  RFl  bias 
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Sqrt(Pb/(Ti+Zr))  =+0.59938+3.76765E-003  *  RF  Power-0.075642*  RFl  bias 


The  above  equations  were  then  used  to  predict  the  thickness  at  various  conditions  (shown  in  red  box  below) 
and  there  is  good  match  between  the  actual  and  predicted  values  for  film  thickness.  Refractive  index  depends 
only  on  the  oxygen  content  of  the  film  and  absorption  coefficient  is  affected  by  any  plasma  variation  (Table  3). 
But  these  ANOVA  models  are  still  useful  in  predicting  the  conditions  to  get  high  deposition  rate  PZT  thin 
films  without  sacrificing  its  optical  properties. 


Table  3:  Validation  (top  box)  and  Prediction  (bottom  box)  of  optimal  operating  conditions  for  a  RF  sputtering 
of  PZT  thin  films. 


U£in£  these  settings  in  equations  on  previous  slide  to  predict 


Thickness 

Refractive  ml 

apx  (nj 

ABstJTT 

Actual 

Predicted 

Actual  Predicted 

Actual 

Predicted 

1  3mT/l2Ar/a00w 

Pt 

1410.39 

1345.19 

2.2758 

2.083 

8.12105 

3.40E.05 

52/48 

^  3mT/l2Ar/200W  +  10WRFl 

•Pt 

708.62 

688,7 

2.0853 

2.083 

1,64E‘12 

1.46E41 

S2/48 

L  3mT/U  +  isccms  02 

i 

Pt 

918.31 

856.SS 

2.028S 

1.89 

S.S7E-03 

0.007673 

A  2“^  DOE  was  conducted  to  ensure  that  we  have  more  robust  models  around  the  preferred  operating  point 
and  the  resultant  responses  are  shown  in  Table  4. 

Table  4:  2"^^  full  factorial  statistical  design  of  experiment  in  3  process  variables  and  the  measured  VASE 
responses. 


Sid 

Run 

Stock 

Factor  1 

A:  02  flow 

seem 

Factor  2 
B:RF  Power 
W 

Factor  3 
C:RF1  bias 
W 

Response  1 
tf^cicness 

A 

Response  2 
n@633ntn 

Response  3 
k@633nm 

3 

1 

Btockl 

000 

200  00 

000 

117971 

2.3469 

0  010492 

7 

2 

0tock1 

0.00 

200.00 

2.00 

1025.72 

2. 2847 

0.0013355 

8 

3 

0tock1 

2.00 

200  00 

2.00 

BOSS 

2.1958 

0.011107 

1 

4 

0tock1 

0.00 

150.00 

0.00 

817.54 

2.2963 

O.OOS3137 

6 

S 

Block  1 

2  00 

150  00 

200 

25621 

3.2466 

0.020467 

S 

6 

0tock1 

1.00 

175.00 

jv  1.00 

556.29 

2.0687 

0.0097767 

2 

7 

Btockl 

200 

150.00 

^  000 

390.71 

2.0197 

0  044156 

4 

8 

Btocki 

2.00 

200.00 

0.00 

719.63 

2.0161 

0.0033381 

5 

9 

Btockl 

000 

150  00 

200 

64292 

2.2615 

0.0020284 

The  results  of  analysis  of  variance  of  VASE  data  gave  us  the  following  models  for  responses  (this  time  only 
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absorption  coefficient  had  to  be  transformed  for  more  linearity  in  the  model).  Noticeably,  most  of  the 
interaction  terms  have  dropped  out  of  the  model.  The  expected  effect  of  plasma  variables  on  refractive  index 
and  absorption  coefficient  is  more  evident  and  realistic. 

Thickness  =  -157.05500-224.18000*  02  flow+6.61790*RF  Power-84.60500  *  RFl  bias 


n@633nm  =  +2.32160-0.15185  *  02  flow-0.024250  *  RFl  bias+0.062950  *  02  flow  *  RFl  bias 

1.0/(k@633nm)  =  +104.38522-21.11945  *  02  flow  +207.25991*  RFl  bias-  101.80498  *  02  flow  *  RFl  bias 

Based  on  the  SDE’s,  we  targeted  a  1  micron  PZT  film  with  Pb/(Zr+Ti)  ratio  of  1  and  low  absorption 
coefficient.  SEM  cross-section  in  Fig.  19  shows  a  film  of  989  nm  whilst  our  VASE  measurement  was  903nm. 
The  slight  discrepancy  is  due  to  low  magnification  of  the  SEM  picture.  On  careful  perusal,  one  can  observe  a 
thin  initiation  layer.  This  is  a  high  deposition  rate  thin  layer  of  80- 1  OOnm  PZT  that  will  provide  nucleation  for 
creating  crack  free  films. 


Jlin  AO-  7min  EhT-SWHV 

I  I  M*  >  1  A  * 


Figure  19:  SEM  cross  section  of  the  PZT film  over  a  260nm  sputtered  Pt/  5  Onm  Ti/500nm  Si02/Si  substrate. 


We  also  conducted  an  EDX  scan  of  the  PZT  thin  film  and  evidently  as  shown  in  Fig.  20,  we  did  obtain  the 
Pb/(Zr+Ti)  ratio  of  1.  The  individual  elemental  concentrations  [Pb],  [Zr]  and  [Ti]  were  incongruent  with  the 
ANOVA  models  as  the  EDS  analysis  equipment  and  software  were  dissimilar.  Earlier  EDX  data  was 
quantitative  whilst  this  is  qualitative  due  to  the  standard  less  nature  of  the  analysis. 


Figure  20:  EDX  stoichiometry  results  of  the  PZT  film  in  the  blue  box. 
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Figure  21  shows  the  ferroelectric  response  of  the  deposited  films  exhibiting  remnant  polarization  of  Pr  =  28.04 
pC/cm^  and  22.38  pC/cm^  which  is  close  to  that  obtained  for  a  sol-gel  seed  layer  textured  RF  sputtered  thin 
film.  The  coercive  field  was  Ec=  444.7  kV/cm  &  475.6  kV/cm. 


Figure  21:  Ferroelectric  hysteresis  loops  for  -0.9pm  thick  PZT  sputtered  film. 

Two  measurements  as  shown  in  Fig.  22  were  attempted  to  isolate  the  effect  of  substrate  clamping  -  in  plane 
measurement  depicted  as  ‘dot2dot’  and  out  of  plane  depicted  as  ‘top2bottom’.  The  dot2dot  loop  is  more 
symmetric  due  to  the  lack  of  residual  stress  and  substrate  clamping  effects. 


Figure  22:  Test  probe  configuration  for  in-plane  vs.  out-of-plane  measurements. 

The  results  of  the  2"^^  DOE  Analysis  of  Variance  of  VASE  data  gave  us  the  following  models  for  responses: 

Thickness  =  -157.05500-224.18000*  02  flow+6.61790*RF  Power-84.60500  *  RFl  bias 

n@633nm  =  +2.32160-0.15185  *  02  flow-0.024250  *  RFl  bias+0.062950  *  02  flow  *  RFl  bias 

1.0/(k(5)633nm)  =  +104.38522-21.11945  *  02  flow  +207.25991*  RFl  bias- 101.80498  *  02  flow  *  RFl  bias 

We  also  conducted  EDX  scans  of  the  PZT  thin  films  generated  by  these  2  SDE’s.  The  ANOVA  models  for  the 
same  are  shown  below.  But  we  decided  to  proceed  further  and  investigate  a  method  to  correlate  the  optical 
properties  of  these  films  with  the  elemental  composition  (as  in  Figure  22). 

[Pb]  =-5.15252+17.71069  *  02  :Ar  ratio+0.092922  *  RF  Power+0.34365*  RFl  bias-0.28374  *  02  :Ar  ratio  *  RF 

Power+3.57606*  02  :Ar  ratio  *  RFl  bias-7.58598E-003*  RF  Power  *  RFl  bias 
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Ln(Ti)  =+1.09901-2.93732  *  02  :Ar  ratio 


[Zr]  =-0.12168+5.48020  *  02  :Ar  ratio-0.13771  *  Pressure+0.025982  *  RF  Power+0.12614  *  RFl  bias-0.10126 
*  02  :Ar  ratio  *  RF  Power+1.12925  *  02  :Ar  ratio  *  RFl  bias-1.91314E-003*  RF  Power  *  RFl  bias 

Sqrt(Pb/(Ti+Zr))  =+0.59938+3.76765E-003  *  RF  Power-0.075642*  RFl  bias 


Figure  22:  Typical  Optical  properties  vs.  wavelength  for  a  particular  PZT  film  and  its  elemental  composition 
as  measured  by  Energy  Dispersive  X-ray  Analysis. 


We  plotted  and  computed  Tauc  and  Wemple-DiDomenico  parameters  from  the  optical  dispersion  data  (Figure 
23).  We  then  regressed  any  correlation  between  elemental  composition  and  these  parameters.  But  we  needed 
another  parameter  to  equate  with  the  number  of  unknowns  in  the  resultant  model  equations.  From  Figure  24, 
electronic  polarizability  affects  the  optical  properties  of  a  thin  film  and  so  we  computed  it  using  3  different 
approaches  (as  plotted  in  Figure  25). 


•  Tauc  plots 
®  Derive 


•  ^Vemple-DiDoniemco  relationship 
•  Derive  Ej.  E, 


Figure  23:  Tauc  and  Wemple-DiDomenico  plots  for  T^  SDE  PZT  thin  films. 
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Figure  24:  Electronic  Polarizability  and  Optical  dielectric  response  of  a  thin  film 


Figure  25:  Electronic  polarizabilities  calculated  using  3  different  approaches. 

Using  the  Tauc,  Wemple-DiDomenico  and  Polarizability  correlations  to  composition,  we  can  predict  the 
composition  of  a  thin  film  from  the  same  system  by  substituting  the  measured  optical  dispersion  parameters. 
The  model  equations  for  the  SDE  and  the  model  predictions  using  this  methodology  are  shown  below.  This 
methodology  translates  into  a  lab  scale  non-destructive  compositional  prediction  model  for  multi  component 
thin  films. 

*  Itsmg  fcjr  Ec  or  Ed  or  Eg  Irnmi  saniplr  and  prcnoiislv  derived 

correUlioii.'i  (sliou  n  Kdott  ),  solve  I'or  .  Ijj  Jnd  t',-j 

^ .  _  fpbf^Pb  ^ 

V  2  (fzr^Zr  +  /rt^>c  +  fo^o) 

ft,  -  -0.132151  +  0.0375312 
fzr  =  -0.066168  +  0.0193297  «  Eo 
fpi,  =  0.3855925  -  0.0586112  -  Eo 
a„  -  -6S9e  -  25  +  4.0  le  -  24  t/ 
a^rr  =  1.82  le  -  24  +  0,452158  ^ 

^  1.348e  -  24  +  07085994 

Figure  26:  The  correlations  (left)  and  the  predictability  (right)  of  the  new  methodology. 

4.3  Magnetic  Field  Sensor  Fabrication 

The  single  layer  transformer  structure  is  shown  in  Fig.  27.  It  consists  of  a  nickel  ferrite  (NFO)  dot  over  a 
patterned  platinum  electrode  which  in  turn  is  over  a  patterned  or  un-pattemed  PZT  thin  film.  A  diffusion 
barrier  layer  might  be  incorporated  between  the  PZT  and  Si  so  as  to  prevent  Si02  formation.  It  could  a  layer  of 
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Pt  which  can  serve  as  the  ground  for  both  input  (ring)  and  output  (dot)  electrodes. 


Urcuuirii 


Figure  27:  Schematic  of  a  single  layer  transformer-based  sensor  structure. 

We  have  utilized  photo  etched  metal  shadow  masks  to  deposit  patterned  thin  films  (see  Table  below  for 
process  sequence).  We  also  used  magnets  (shown  in  blue  in  Figure  28)  to  protect  the  underlying  layer’s 
electrical  pads  from  deposition. 
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Figure  28:  Shadow  mask  processing  using  metal  shadow  mask  and  magnets  (to  protect  electrical  pads  from 
deposition). 


We  had  to  develop  eustom  testing  eapability  for  testing  our  ME  deviees.  We  proeured  a  Agilent  E4991A 
Impedanee  Analyzer,  Jmiero  Teehnology  probe  station  and  GGB  Pieoprobes  with  speeial  wire  tethers  (to 
ereate  eommon  ground  between  input  and  output).  Figure  29  shows  the  measurement  system. 


b)  c) 


*  RFanal\T:er  •  ImicrnTecbioiogy  kRN  -()9S 

*  JmicToTi’rhnology  LMS-2709  PosttiomTs 

Priibc,'  St  jtion  *  GGB  Piroprnbe  ECP I S  GS  ^^,St  1 


Figure  29:  Electrical  connections  and  equipment  for  High  Frequency  ME  testing  -  a)  DUT  (device  under 
test),  b)  impedance  testing  schematic,  c)  Gain  testing  schematic,  d)  test  bench  and  e)  probe  tip  with  special 
wired  tethers  for  common  ground. 


We  first  fabrieated  the  multilayer  piezoeleetrie  transformer  based  on  tape-easting  teehnique  and  low 
temperature  eo-firing  proeess.  The  advantage  of  this  teehnique  is  low  eost  and  availability  for  industry 
manufaeture.  With  optimized  inner  eleetrode  eomposition  and  eo-firing  proeess,  we  were  sueeessful  in 
aehieving  high  density  and  eraek  free  eo-fired  PZNT  (0.8PZT  (52:48)-0.2PZN)  multilayer  transformer  at 
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930°C.  On  the  co-fired  multilayer  transformer,  we  attaehed  two  layers  of  Met-glass  on  both  side  of  the  output 
section  with  dimension  of  4x4  mm^  to  fabricate  the  laminated  ME  transformer  as  shown  in  Fig.  30.  With  Met- 
glass  attached  on  the  output  area,  we  were  able  to  tune  the  performance  of  the  transformer  under  external  DC 
magnetic  bias  as  shown  in  Fig.  31  and  Fig.  32. 


Output 


Figure  30:  Schematic  of  ME  transformer  structure. 


Figure  30:  (a)  ME  transformer  voltage  gain  as  a  function  of  frequency  at  varying  external  DC  magnetic  field 
from  0  Oe  to  3000  Oe,  (b)  Variation  of  resonance  frequency  as  a  function  of  external  DC  magnetic  field,  (c) 
Change  in  maximum  voltage  gain  with  respect  to  change  in  the  DC  magnetic  field. 


Votage  gain  without  external  magnetic  field 

Figure  31:  ME  transformer  characterization  showing  the  change  of  voltage  gain  in  the  frequency  range  of 
110  kHz-113  kHz  with  and  without  applied  DC  magnetic  field. 
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When  the  transformer  was  operated  under  external  magnetie  field  range  of  400  Oe  to  3000  Oe,  the 
value  of  maximum  voltage  gain  and  its  working  frequeney  were  found  to  inerease  first  with  inereasing  Hdc  and 
then  deerease  with  inereasing  Hdc-  Also,  we  ean  notiee  that  the  voltage  gain  tunability  (Gh-Gq/Go)  of  the 
transformer  is  small  at  low  magnetie  field  below  1000  Oe  but  inereased  at  high  magnetie  field.  The  trend  of 
resonanee  frequeney  shift  was  found  to  be  as  expeeted  from  the  equivalent  eireuit  models.  The  data  elearly 
demonstrates  our  hypothesis  of  tunable  performanee  of  ME  transformer  under  external  magnetie  field. 

Graded  transformer:  Sinee  the  piezoeleetrie  transformer  ean  be  seen  as  a  eombination  of  aetuator  (input) 
and  transdueer  (output),  we  experimented  with  a  new  transformer  strueture  with  different  materials  at  input 
and  output  seetion.  Materials  whieh  possess  higher  dss  are  suitable  for  aetuator  and  those  with  higher  Qm  are 
suitable  as  a  transdueer.  This  design  provides  the  potential  of  aehieving  better  performanee  and  new  strueture 
of  ME  tunable  transformer.  Figure  32  shows  the  strueture  of  graded  transformer.  This  figure  shows 
eomparison  of  voltage  gain  and  effieieney  performanee  between  single  phase  piezoeleetrie  transformer  and 
graded  transformer.  Further  analysis  has  to  be  eondueted  on  this  new  design  with  aim  of  improving 
piezoeleetrie  transformer  performanee. 


Figure  32:  Graded  piezoelectric  transformer  with  ring-dot  structure. 


4.4  Hybrid  deposition  process 

We  developed  hybrid  deposition  proeess  to  aehieve  high  quality  films  by  eombining  pulsed  laser 
deposition  (FED)  and  sol-gel.  Pulsed  laser  deposition  has  high  deposition  rate  and  flexibility  in  target  material 
but  as-deposited  films  are  not  well  textured.  So  by  taking  the  advantage  of  orientation  eontrol  in  sol-gel 
proeess  and  high  deposition  rate  of  FED  proeess,  we  expeet  to  aehieve  high  quality  textured  film.  The  sol 
preparation,  gelation  proeess  and  perovskite  erystallization  proeess  were  optimized  for  (100)  orientation. 

Pulsed  laser  deposition  with  parameter  (9E-4  Torr  ambient  oxygen  pressure,  90  min,  room  temperature)  was 
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conducted.  Film  were  grown  at  room  temperature  on  as  deposited  PZT  seed  layer  (sol-gel).  In-situ  annealing 
temperature  was  varied  in  the  range  of  250°C  to  650°C.  Crystallization  was  determined  by  using  X-ray 
diffraction  (XRD)  as  shown  in  Fig.  33.  We  can  notice  that  this  process  provides  <100>  preferred  orientation. 
The  fraction  of  texturing  increased  20%  from  58.7%  (seed  layer)  to  79.12%  (2“^  deposited  film)  by  using  this 
hybrid  process. 

PZT  (100)  Sol-gel  &  PLD  hybrid  deposition  Fraction  of  (100)  texturing 


79,12%  (optimized) 


Figure  33:  XRD  pattern  of  <100>  textured  PZT  film  on  platinized  Si  via  hybrid  deposition. 


Using  the  films  deposited  through  this  combinatory  process,  we  started  to  fabricate  the  unipoled 
transformer  structure.  Figure  34  shows  a  series  of  image  of  the  first  prototype.  Ring-dot  silver  electrode  was 
printed  by  aerosol  jet  deposition  system.  Further  investigation  on  voltage  gain  will  be  conducted  in  the  near 
future. 


Figure  34:  Unipoled  PZT  thin  film  transformer 
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Figure  35  shows  the  developed  hybrid  deposition  proeess  by  eombining  pulsed  laser  deposition  (PLD)  and 
sol-gel.  High  deposition  rate  ean  be  aehieved  by  this  proeess  for  thiek  film  growth.  In  addition,  high  fraetion 
of  texture  ean  be  obtained  via  optimized  deposition  parameter. 


»  »  M  H  TO  IP  H  4  n  «  TP  n  »  i«  iP  PP 


2*  funs? 


Figure  35:  Schematic  of  hybrid  deposition  process,  XRD  evolution  of  hybrid  deposited  PZT  films  under 
varied  annealing  temperature,  microstructure  and  piezoelectric  properties  of  (100)  highly  textured 
Pb(Zr0.6Ti0.4)03  thick  film. 


Hetero-structure  growth  ofBaTiOfBTO)  &  CoFeyOfCFO) 

We  studied  the  growth  and  mierostrueture  of  BaTiOs  (BTO)  thiek  films  on  platinized  silieon  substrate 
synthesized  by  pulsed  laser  deposition  as  shown  in  Fig.  36.  We  investigated  the  evolution  of  mierostrueture 
and  piezoeleetrie  response  for  as-grown  BTO  film  as  a  funetion  of  thiekness.  Interestingly,  the  surfaee 
morphology  of  the  films  adopted  the  symmetry  of  underlying  layer  and  evolved  from  pyramid  to  hexagram 
shape  eolumnar  strueture  with  (lll)-preferred  orientation  as  the  film  thiekness  was  inereased.  The 
piezoeleetrie  properties  were  also  found  to  inerease  with  thiekness.  We  expeet  that  the  strueture-property 
relationship  deseribed  for  the  thiek  BTO  piezoeleetrie  film  in  this  study  will  strengthen  the  integration  of  BTO 
with  silieon. 

For  thin  film  ME  tunable  transformer,  it  is  advantageous  to  study  the  magnetoeleetrie  heterostrueture 
that  offers  the  potential  of  magnetie  or  eleetrie  field  indueed  tunabality.  We  have  not  only  developed 
eonventional  2-2  type  strueture  to  study  the  interaetion  between  two  different  phases,  but  also  ereated  a  new 
magnetoeleetrie  strueture  whieh  might  possesses  better  performanee  on  ME  eoupling.  Figure  37  shows  TEM 
image  of  ion  bombardment  assisted  as-grown  BaTi03/CoFe204  heterostrueture  and  its  ferroeleetrie  properties. 
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Ion  bombardment  induced  defect  were  found  as  a  result  of  merged  columnar  structure  which  could  provide  us 
a  new  method  to  tune  ferroelectric  properties  via  modulation  of  defect  density. 


Figure  36:  Morphology  evolution  and  piezoelectric  properties  of  as-grown  BaTiOs  thick  film. 


The  interaction  between  these  two  phases  was  studied  with  varying  deposition  temperature  and  BTO 
thickness.  Deposition  parameter  for  single  layer  BTO  and  CFO  film  was  optimized  as  lOOmTorr  ambient 
oxygen  pressure  and  deposition  temperature  of  700-850°C.  Further  study  on  effect  of  variables  such  as  BTO 
thickness  and  growth  temperature  was  conducted  as  shown  in  Fig.  38.  Grain  size  and  thickness  were  found  to 
be  quite  sensitive  to  the  BTO  thickness  and  growth  temperature.  Triangular  topology  was  found  to  dominate 
grain  shape  SEM  which  reflect  film  growth  on  (111)  oriented  Pt  as  shown  in  Fig.  39. 
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Figure  39:  SEM  image  of  IBAD-CFO(left)  and  BTO  on  as-grown  IBAD-CFO  film. 


We  have  quantified  magnetization  and  polarization  as  a  function  of  ion  bombardment  time  as  shown 
in  Fig.  40.  As  expected,  saturation  magnetization  decreased  with  increasing  ion  bombardment  time.  The 
reason  for  this  trend  was  related  to  decrease  in  CFO  film  thickness  with  increase  in  ion  bombard  time.  P-E 
loop  in  this  heterostructure  shows  better  hystersis  than  simply  BTO/CFO  and  BTO  films.  The  reason  for  this 
change  still  under  investigation.  We  also  measured  ferroelectric  properties  change  as  a  function  of  magnetic 
field  as  shown  in  Fig.  41.  We  can  clearly  notice  the  variation  of  ramanent  polarization  change  which 
demonstrates  the  existence  of  coupling  between  magnetic  and  ferroelectric  components. 
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Figure  40:  (a)  M-H  loop  of  IBAD-CFO  films;  (b)  P-E  loop  of  BTO/IBAD-CFO  films;  (c)  Comparision  of 
magnetization  of  IBAD-CFO  film  and  (d)  Comparision  of ferroelectricity  as  a  function  of  Ion-bombardment 
time. 
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Figure  41:  Ramanent  polarization  change  as  a  function  of  magnetic  field  for  as-grown  IBAD  heterostructures. 


Complex  3D  ME  composites 

In  order  to  break  the  conventional  phase  connectivity,  we  investigated  novel  composite  structures  that  can 
overcome  the  limitations  imposed  by  material  symmetry  as  shown  in  Fig.  42.  In  this  process,  we  start  with  a 
pattern  deposited  by  aerosol  deposition  (AD)  on  a  platinized  silicon  substrate  and  then  use  pulsed  laser 
deposition  (PLD)  to  grow  multilayers  of  BaTiOs  and  CoFe204.  After  in-situ  annealing  at  high  temperatures, 
the  deposited  films  self-assemble  into  complex  structures.  Lastly,  we  use  sol-gel  to  achieve  homogeneous  top 
surface.  A  combination  of  these  three  techniques  can  be  used  to  fabricate  large-area  arrays  of  magnetoelectric 
structures.  This  synthesis  method  not  only  helps  in  achieving  a  well-ordered  arrangement  of  the  nanoscale 
magnetoelectric  structures  but  also  simplifies  the  classical  top-down  patterning  and  multi-step  processing.  We 
have  been  focusing  on  synthesis  of  these  structures  and  further  would  like  to  conduct  detailed  investigation  on 
understanding  of  their  static  and  dynamic  ferroelectric  and  ferromagnetic  behavior  at  nanometer  scale. 


ferroelastic  domain  structures,  and  magnetoelectric  coupling. 
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Figure  42:  Design  and  synthesis  of  complex  3D  ME  composites. 
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Using  this  process,  we  designed  and  synthesized  a  new  structure  of  CFO/PZT  ME  composite  array  with  multi¬ 
orientation  as  shown  in  Fig.  43.  With  highly  textured  buffered  array  pattern  (aerosol  jet  deposited  PZT)  and 
optimized  growth  condition  of  CFO/PZT  heterostructure,  we  successfully  obtained  composite  array  with 
multiorientation.  CFO  inside  the  pattern  shows  highly  (111)  prefered  orientation  while  outside  the  pattern 
shows  hihgly  (100)  prefered  orienteation  with  smaller  grain  size. 
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Figure  43:  CFO/  PZT  ME  composite  array  with  multi-orientation. 


Interestingly,  MFM  image  of  as  grown  film  shows  different  magnetic  response  closely  related  with 
grain  size  and  orientation  at  different  area  of  the  array.  CFO  with  larger  grains  possesses  larger  magnetic 
response,  while  small  grains  possess  relative  smaller  response.  The  reason  for  the  response  variation  might  be 
grain  size  and  orientation  difference.  This  provide  us  the  information  that  we  might  be  able  to  tune  magnetic 
response  by  modifying  grain  size  and  orientation  during  heterogeneous  growth.  Moreover,  we  can  obverse 
higher  piezoresponse  on  small  CFO  grains  which  show  relative  small  magnetic  response  inside  the  pattern. 
This  might  be  explained  through  magnetoelectric  coupling.  With  applied  bias  voltage,  PZT  film  below  will 
transfer  strain  to  the  magnetic  layer  which  in  return  shows  high  piezoreponse  amplitude  on  the  surface. 


Figure  44:  Magnetic  force  microscopy  (MFM)  images  and  Piezoresponse  force  microscopy  (PFM)  images  of 
as  grown  film. 

Conclusions 

We  have  developed  textured  PZT  thin  films  using  two  different  deposition  techniques,  sol-gel  and  RF 
sputtering.  We  have  characterized  these  films  for  their  composition,  structure  and  crystallinity.  An  upshot  of 
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the  detailed  characterization  was  the  development  of  prediction  models  for  texturing  of  PZT  sol-gel  thin  films, 
an  understanding  of  the  analytical  techniques  that  can  discriminate  between  the  highly  textured  films, 
optimization  of  the  RF  sputtered  PZT  thin  film  properties  using  design  of  experiments  methodology  and 
finally  the  establishment  of  a  lab  scale  non-destructive  compositional  analysis  methodology  for  PZT  RF 
sputtered  thin  films  using  Ellipsometry.  We  have  conceived  a  unipoled  thin  film  transformer  device  utilizing 
the  existing  toolset  at  VT’s  MicrON  clean  room  facility  (for  mask  alignment)  and  our  sputtering  capability. 
We  have  fabricated  single  layer  transformer  structure  and  characterized  its  performance. 
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